New Luttinger liquid physics from photoemission on Lio.gMogOiy 



O 

o 

(N 

< 

(N 






I 

o 
o 



> 

m 
o 
in 

o 

o 

-(— • 



O 

o 



X 



Feng Wang,i J. V. Alvarez,2.E|s.-K. Mo,i J. W. Allen,i G.-H. Gweon,^ J. He,"^ R. Jin,^ D. Mandrus,-*- ^ and H. Hochst^ 

'Randall Laboratory of Physics, University of Michigan, Ann Arbor, Michigan 48109, USA 

^Departamento de Fisica de la Materia Condensada, 

Universidad Autonoma de Madrid, 28049 Madrid, Spain 

^Department of Physics, University of California, Berkeley, California 94720, USA 

^'Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37996, USA 

^ Solid State Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 

^Synchrotron Radiation Center, University of Wisconsin, Stoughton, Wisconsin 53589, USA 

(Dated: February 6, 2008) 

Temperature dependent high resolution photoemission spectra of quasi-1 dimensional Lio.gMoeOi? 
evince a strong renormalization of its Luttinger liquid density-of-states anomalous exponent. We 
trace this new effect to interacting charge neutral critical modes that emerge naturally from the 
two-band nature of the material. Lio.gMoeOir is shown thereby to be a paradigm material that is 
capable of revealing new Luttinger physics. 

PACS numbers: 71.10.Pm, 71.10.Hf, 79.60.-i 



Condensed matter physics owes much to research on 
paradigm materials, i.e. materials that typify a basic 
phenomenon and often a whole class of materials. The 
best paradigm materials not only provide a proof of ex- 
istence in nature, but also reveal new physics of the phe- 
nomenon. E.g., research on silicon has shaped the under- 
standing of covalently bonded semiconductors. This pa- 
per reports research on Lig gMogOiy that not only solid- 
ifies its status as a paradigm quasi-1 dimensional (QID) 
material showing Luttinger liquid (LL) physics, but goes 
beyond previous work on this or any other QID material 
to reveal new LL physics as well. 

The LL theory |]| provides an understanding of the 
ID orthogonality catastrophe, in which single-electron 
behavior cannot be adiabatically continued from a ID 
system of free electrons when electron-electron interac- 
tions are included. This striking phenomenon is di- 
rectly observable in angle integrated photoemission ex- 
periments 131 as a nearly complete power law suppres- 
sion of the single particle density of states (DOS) at the 
Fermi energy Ep. LL theory also characterizes (1) the 
low-energy excitation spectrum as consisting entirely of 
independent, linearly dispersing spin and charge density 
fluctuations (spinous and holons, respectively) and (2) 
the various correlation functions as displaying power law 
behaviors with anomalous exponents, e.g. a for the single 
particle DOS near Ep. For this theory to apply in a real 
solid material, it is essential to have (i) high anisotropy 
and (ii) linearity of the dispersion of the underlying bands 
over a substantial energy range away from Ep. The im- 
portance of (ii), e.g. in setting the energy scale for ob- 
serving LL behavior, tends to be overlooked. 

Among all the few QID solids |3| for which LL behav- 
iors are claimed, Lio.gMogOiy is unique as a metal with 
a well understood band structure i^j that clearly satisfies 
both criteria. Strong anisotropy has been confirmed by 
various experiments ja, |y| and linear dispersion over a 



range w 0.13 eV from Ep is observed in an gle resolved 
photoemission spectroscopy (ARPES) [ili.l9tllollTlllli|. 

which directly measures the single particle spectral den- 
sity as a function of momentum k and energy to. 

ARPES spectra 0, 0, iJH El E3 have shown ex- 
pected holon and spinon 'll'| features near Ep, and also 
a power law DOS when angle integrated. These findings 
are strongly supported qualitatively by recent tunneling 
spectroscopy "iS' that also shows a power law DOS near 
Ep. However, the values of a, « 0.6 and « 0.9, extracted 
from tunneling and ARPES, respectively, are quite dif- 
ferent. The only hope of reconciling this difference lies 
in the differing T ranges, 5-55 K, and 250-300 K of the 
two experiments, respectively. Unfortunately, as we show 
below, such a large T-renormalization of a is inconsis- 
tent with Lio.gMogOiv transport properties implying T- 
independent a within usual one-band LL theory. In this 
Letter we report new ARPES directly showing that the 
large T-renormalization of a does indeed occur and we 
show that its origin lies in interacting charge neutral crit- 
ical modes that emerge naturally for Lio.gMogOiy be- 
cause two bands cross Ep. The theory also predicts the 
possibility of observing two distinct kinds of a(T) be- 
havior on various sample surfaces, as may actually have 
been observed. This is both the first observation and 
the first theory of such a strong anomalous exponent T- 
dependence. 

ARPES experiments over a wide T-range 15-300 K, 
with photon energies 24 eV and 30 eV, were performed 
at the Ul-NIM undulator beamline of the Synchrotron 
Radiation Center of the University of Wisconsin. Sin- 
gle crystal samples were grown by the temperature gra- 
dient flux method and were cleaved in-situ at pressure 
< 1 X 10~i° torr for measurement. A Scienta SES2002 
spectrometer with an acceptance angle of 14° was used 
for all data. Freshly evaporated gold was measured be- 
fore and after the experiment at temperature ^ 15 K, 
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FIG. 1: ARPES spectra of Lio.gMoeOir measured at T = 200 
K, hu — 30 eV. (a), (b): Energy versus momentum spectra 
taken along T ~Y and off- normal direction, respectively, (c): 
A schematic view of the 1-D FS with the k space paths of (a) 
and (b) indicated, (d): Stack view of (b), with dashed lines 
highlighting the band dispersions. 

for a careful calibration of the Fermi enegy as well as the 
instrumental energy resolution, AE « 18 meV FWHM. 
We have collected extensive data for samples cleaved at 
different temperatures and with the measurement tem- 
perature both increasing and decreasing. 

The ARPES spectra in Figs. 1(a) and 1(b) were taken 
with angular resolution « 0.3° and the analyzer slit per- 
pendicular to the QID Fermi surface (FS) to show the 
band dispersions. The associated k-space paths are il- 
lustrated in Fig. 1(c). Lio.gMogOiy has 4 Mo 4d bands 
near Ef, labeled A - D in Fig. 1(d). Two of them (C 
and D) become degenerate at binding energy ss 0.13 eV, 
then disperse linearly [ij and cross Ep together to de- 
fine the QID FS. Along path (a), the intensity of band 
D is very suppressed by the ARPES matrix element. 
However, both C and D are clearly seen along path (b). 
The measured band structure is in good agreement with 
that found previously by ARPES and qualitatively with 
the band structure calculations. The same spectra were 
taken at various temperatures down to 15 K and show 
no apparent change of the bands. 

Angle integrated spectra were taken for the same ge- 
ometries as for Figs. 1(a) and 1(b). Fig. 2(a) shows 
a set of spectra for selected temperatures, with the in- 
set showing a wide range spectrum for 30 K. All spectra 
are normalized to the intensity between binding energy 
0.2-0.8 eV, which is mainly from bands A and B in Fig. 
1. At all T the power law suppression of the near E-p 
DOS is evident. As T decreases there is much sharpen- 
ing, a great reduction of the tail above Ep and a strong 
increase of the curvature. This sharpening is only partly 
a reduction of thermal broadening. There is also a signif- 
icant decrease of the power law exponent a as we have 



FIG. 2: (a): T-dependent spectra angle integrated along path 
for Fig. 1(a). Only selected temperatures are shown for clar- 
ity. Inset shows the wide energy range spectrum at 30K. (b), 
(c), (d): Fitting with finite T LL spectral function. Inset in 
(b) shows the details of the fitting between 0.04-0.07 eV. 

confirmed and quantified by comparing the experimental 
spectra with the LL DOS using the finite temperature 
theory of Ref. [ij]. The calculated theory spectra were 
convolved with a Gaussian with the width set by the 
instrumental energy resolution AE, and were scaled as 
needed because the relative experimental intensities for 
various T are not well calibrated. As illustrated in Figs. 
2(b)-(d) the experimental spectral shape at each T can be 
very well fit by the theory up to binding energy 0.105 eV, 
above which the A band weight starts to come in. The 
a value is T-dependent as plotted in Fig. 3(a). Consis- 
tent with previous photoemission and tunneling, a starts 
from « 0.9 at room temperature, decreases sharply be- 
tween 300 K and 150 K, then slows down and saturates 
at ^ 0.6 between 150 — 50 K. At even lower temperature, 
there is a slight revival of a, as documented in Fig. 2,(b) 
and (c). Over the whole range, a is never less than 0.6. 
Fig. 3(a) also shows a{T) derived from spectra taken 
with the geometry of Fig. 1(b). 

Our interpretation of this q:(T) behavior, which yields 
the solid theory curve of Fig. 3(a), flows from the ex- 
cellent agreement, found here and previously [ll|, Il2l |. 
between finite-T LL [lj| spectral theory and measured 
photoemission spectra in the full range of T over which 
a (which itself is very large) renormalizes. Therefore, the 
most suitable starting point for our calculation is an effec- 
tive theory based on the LL where T establishes the low- 
energy cutoff. However, in a single-band spin rotational 
invariant LL any variation of a comes from the renor- 
malization of Kp, the Luttinger parameter in the charged 
sector through the relation a = {Kp + K~^ — 2)/4. The 
driving force of that renormalization may be (i) Umk- 
lapp scattering, present only when Akp = G where kp 
is the Fermi vector and G is a reciprocal lattice vector 



or (ii) for an incommensurate system, the interplay of 
weak disorder and backward scattering, which recouples 
charge and spin [l^. Case (i) is forbidden here because 
kp in Lig gMogOiy is incommensurate |^. Further, in 
both cases any temperature change in a is necessarily 
caused by renormalization in the charged sector and so 
must be correlated with strong transport anomalies in 
exactly the same range of temperatures. In contrast, we 
have found that the resistivity for Lio.gMoeOiy displ ays a 
nearly perfect power-law in T for almost one decade |lfi| . 
This lack of scaling relation between the spectral and 
transport exponents suggests that the renormalization 
process is taking place in a neutral sector. Both our pho- 
toemission experiments and band theory calculations |j| 
show two bands crossing the Fermi level, which offers 
a natural origin for such neutral modes as the fluctua- 
tions in the difference of density between the two bands 

P-i^) = {pix)^ + p{x)ci) - {p{'x)DT+p{x)Di)- Follow- 
ing the standard bosonization process we can decouple 
the fixed-point Hamiltonian in the form: 
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where Il{x)af3 and (f>{x)ai3 are the conjugate bosonic 
fields, Vaf3 the sound velocities and Kaf3 the Luttinger 
parameters(/3 = +/— superposition/difference, a — p/a 
charge/spin density). The p_ and cr_ modes can be pic- 
tured as the ID-allowed bosonic excitations into which 
an interband exciton can decay. Having two extra modes 
with different velocities, we expect wide ARPES line- 
shapes as has actually been observed in Lio.gMogOiy [llj. 



and now a 



(A'p+ + K;^ + Kp^ + K^l - 4)/8. 
Let us consider now how these critical modes can in- 
teract via marginal interactions (those that can renor- 
malize critical exponents). The interactions compatible 
with the symmetries represent intra and inter-band spin- 
backscattering and pair-tunneling processes between the 
two bands. A convenient form in which to write the 
renormalization group equations is: 

ga-Vi + {d- ga+}V2 , 
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(ff'T+) 
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g^j. represent spin-backscattering process, Vi interband 

pair-tunneling and d = 1 — -jy 3ct+- 

Since we have restricted ourselves to spin rotational 
invariant flows, Kcr+ and K^-- do not show up explic- 
itly in the equations. The total charge mode {p+) con- 
tributes to the total value of a but the absence of Umk- 
lapp scattering keeps the charged sector independent and 
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FIG. 3: (a): a(T) extracted from fitting of the spectra in 
Fig. 2(a) (circles) for a complete T series and the same (tri- 
angles) measured for the k-path of Fig. 1(b) on the same 
sample. Typical error bar from fitting « ±0.03. The line is a 
theory calculation with our model (see text), (b): Kosterlitz- 
Thouless flow for V2 = 0, g^- = 0; finite values give the same 
basic structure, (c): q:(T) as Vi increases, showing two differ- 
ent behaviors, (d): q(T) measured on 6 samples, with error 
bars (not shown) similar to (a). Shading highlights two differ- 
ent behaviors much like those in theory curves of panel (c)(see 
text). 

free of significant renormalizations, consistent with the 
T-independence of the resistivity exponent in this range 
of T, and with the low T saturation value of q;(T) being 
greater than 0.6. The lack of commensurability has been 
confirmed by photoemission and the two band structure 
calculations available and reinforces the essential role of 
the neutral mode for a(T). The interchain 3D perturba- 
tions do not seem to play any important role at these 
temperatures, the only evidence of the 3D phy sics being 
the superconducting transition at 1.8 K |5lll6lll7l|. 

We are not interested here in the correlations dom- 
inating at T = which have already been studied in 
detail (see for example Ref. [13] and references therein). 
Rather we focus on the crossover structure at higher tem- 
peratures which in Lio.gMogOiv has great physical signif- 
icance due to the large energy range over which the bands 
crossing Ep have linear dispersion, as pointed out above. 
Fig. Ofb) shows the result of integrating the RG equa- 
tions for varying Vi(0). For V2(0) = (/o— (0) = and 
gcr+ > 0, Eqs. (^ and ^ display literally a Kosterlitz- 
Thoulcss ^9y (KT) flow in the charge density difference 
(p— ), with decoupled gapless excitations for the total 
spin {a+). This behavior is not isolated. To be sta- 
ble in the range of temperature of interest, finite values 
of V2(0) and positive values < ^^-(O) < ^^+(0), it is 

sufficient that A 



d+ a/.9ct+ + 9a- < 0- I3esides captur- 
ing the LL behavior and the instabilities of the two-band 
system, this KT fiow supports simultaneously the correct 



magnetic (Pauli susceptibility |5|,|g]) and spectral (T de- 
pendence of a) phenomenologies. Fig. |3c) shows the 
a(T) corresponding to the flows of Fig. Olb). The the- 
ory curve of Fig. ^a.) gives a good general description of 
the data, thus placing Lio.gMogOir in a region of param- 
eters where the RG flow and the a(T) can be understood 
in terms of the KT behavior, i) The downward renormal- 
ization of a comes from a reduction of the 1/Kp- and a 
simultaneous screening of the V couplings, ii) At small 
values of V1.2 and specially, close to the separatrix ending 
in Kp^ ^ 1, the change in 1/Kp- is very slow (see Eq. 
^). As a consequence, a displays a saturation regime 
over a large range of intermediate temperatures, iii) The 
incipient revival of a takes place when Kp- = 1 and at 
lower temperatures the V's flow to strong coupling. The 
revival takes place at the lowest temperatures we could 
measure so we cannot determine experimentally if it cul- 
minates in a gap opening in the p— sector. If so, that 
crossover and the eventual gap opening might help to un- 
derstand the broad anomaly in the specific heat [a, llg 
that is present in the same range of T. 

The proximity of Lio.gMogOiy to a separatrix suggests 
the possibility that both behaviors of a(T) seen in Fig. 
OJc) might be observed in experiment. As 1/Kp^ de- 
creases the electronic compressibility increases. There- 
fore the p- mode tends to favor phase separation, per- 
haps nucleated by defects in different parts of the sam- 
ple. In the first "phase" the reduction of stiffness is large 
enough to support a revival of a. On the other side of the 
separatrix, the change of a is slower and monotonic. In 
fact, both behaviors may have been observed. Fig. Ol^d) 
plots the a(T) measured for 6 different sample surfaces. 
Each symbol corresponds to one sample except for the 
data repeated from Fig. 3(a). The solid symbols repre- 
sent situations where the sample was cleaved at high T 
and T was then decreased, and the open symbols repre- 
sents situations where the sample was cleaved at 50 K 
and T was then increased after quickly visiting the low- 
est measured T. Despite slight differences from sample 
to sample, our main finding of a decrease and flatten- 
ing of a with decreasing T is observed robustly in all 6 
samples. However, as illustrated by the shadings in the 
figure, it appears that the samples can be separated into 
two groups according to the rapidity of the decrease of 
a(T). The a(T) of one group decreases faster and flattens 
at a higher temperature than occurs for the other group. 
Moreover, the first group tends to show the slight revival 
of ck(T) at low temperature while the slower group does 
not, generally consistent with Fig. OJc). Any variations 
seen in individual data sets beyond these general features 
are deemed to lie within the typical error bars (not shown 
but like those in Fig. El^a)). Further work will be needed 
to fully characterize the two behaviors. 

In summary, T-dependent photoemission of 
Lio.gMogOiv shows that the anomalous exponent a 



has a large renormalization with decreasing T, nicely 
connecting the high and low T values found previously 
in photoemission and tunneling, respectively. New 
theory reveals the essential role of neutral charge 
modes present specifically because there are two bands 
crossing Ep, the absence of which would make the 
observed a(T) impossible to rationalize with a resistivity 
T-dependence that shows no exponent renormalization. 
The possibility of two distinct a(T) behaviors, as may 
have been seen experimentally, is also intrinsic in the 
theory. This is both the first observation of such a 
strong T-renormalization of a(T) and the first theory 
of this unexpected behavior. Moreover, the case that 
Lio.gMoeOiT is the long sought paradigm of LL behavior 
in a solid material has thereby been brought to a new 
level. 
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